Quasinormal modes of black holes in anti-de Sitter space: a numerical study of the eikonal limit 
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Using series solutions and time-domain evolutions, we probe the eikonal limit of the gravitational and scalar- 
field quasinormal modes of large black holes and black branes in anti-de Sitter backgrounds. These results are 
particularly relevant for the AdS/CFT correspondence, since the eikonal regime is characterized by the existence 
of long-lived modes which (presumably) dominate the decay timescale of the perturbations. We confirm all the 
main qualitative features of these slowly-damped modes as predicted by Festuccia and Liu 13] for the scalar- 
field (tensor-type gravitational) fluctuations. However, quantitatively we find dimensional-dependent correction 
factors. We also investigate the dependence of the QNM frequencies on the horizon radius of the black hole 
(brane) and the angular momentum (wavenumber) of vector- and scalar-type gravitational perturbations. 



(N 



CD 

> 

O 
O 

o 

O 



X 
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I. INTRODUCTION 



The characteristic modes of black holes and black branes 
are eigenmodes of these systems which convey important in- 
formation about the background geometry. These systems 
are intrinsically dissipative (energy flows to the event horizon 
and/or towards the spatial infinity) and therefore their eigen- 
modes are not stationary. One calls the eigenmodes of these 
dissipative systems quasinormal modes (QNMs) IH Hi 01 ■ 
Within the anti-de Sitter/Conformal Field theory (AdS/CFT) 
duaUty S i, I?]], these modes serve as an important tool for 
determining the near-equilibrium properties of strongly cou- 
pled quantum field theories, in particular their transport co- 
efficients such as viscosity, conductivity and diffusion con- 
stants m Ht]. Recently, QNMs have also been used to study 
properties like the "meson melting" in D3/D7-brane models 
m [13, mtl and the spectrum of collective excitations of holo- 
graphic superconductors lfl2l [Tsll . 

It is thus of interest to understand the spectrum of black 
holes and black branes in asymptotically anti-de Sitter (AdS) 
spacetimes. Recent studies by Festuccia and Liu (henceforth 
FL) have shed light on this issue, predicting the existence of 
long-lived modes in asymptotically AdS black hole (BH) ge- 
ometries ifTill . Their work shows that the eikonal limit in AdS 
depends sensitively on the relative size of the black hole. For 
small black holes, they find exponentially long-lived modes, 
which can be thought of as modes trapped inside the potential 
barrier If we write oj = cjj^ — iujj for the typical "energy" 
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eigenvalue, these modes take the Bohr-Sommerfeld form 
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/ is the black-hole horizon function (see Eq. ( 12.11 ) below), and 
I is the angular momentum of the perturbation. Accordingly, 
their lifetime is dictated by a "tunneling probability" of the 
form uji cx cxp (— 2 J^'' Qdr), where rf, and Tc < Vb are 
two real zeros (turning points) of Q, and the proportionality 
coefficient is given in II 1511 . This prediction is supported by 
numerical studies for small AdS black holes ifTsll . 

On the other hand, the nature of the long-lived modes for 
large black holes is completely different, since no trapped 
modes are allowed in the large black-hole regime. A con- 
sequence of this fact is that the "Breit-Wigner resonance 
method" used to investigate small AdS black holes fails to 
give reUable results in the large black-hole regime ifTsll . These 
modes are also expected to be long-lived j3 (see Section HTCl 
below). Thus these modes will presumably dominate the BH's 
response to arbitrary perturbations, hence the thermalization 
timescale in the dual CFT. Since their existence may be very 
relevant for the AdS/CFT conjecture, we decided to investi- 
gate numerically these long-lived modes. We use two, well 
established methods to study the QNMs of black holes. One 
is the series solution expansion iflg , [TtI [TsIi . which outputs 
the characteristic frequencies directly and is especially well 
suited for large black holes. By using time-domain methods, 
in particular the scattering of Gaussian wavepackets lfl9l I20I1 
(see also ll2lll ). we confirm and expand the series solution re- 
sults and furthermore show that these modes are excited in 
physically interesting situations. As a by-product, we confirm 
once more that there are no late-time tails in this geometry. 
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II. FORMULATION OF THE PROBLEM 

A. The background spacetime 

BH's in asymptotically AdS spacetimes form a class of so- 
lutions which is interesting from a theoretical point of view 
and central for the study of strongly coupled field theories 
at finite temperature in the gauge/gravity duality framework. 
We are interested in a simple class of non-rotating, uncharged 
d-dimensional Schwarzschild-AdS black holes with line ele- 
ment 



-fdt^ + f-^dr^ 



r'dnl 
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(2.1) 



where f{r) = 1 + r'^/R'^ 



„d-3 



/r 



'^^^ and dfll_2 is the 
metric of the unitary {d — 2)-sphere. The AdS curvature ra- 
dius R is related to the cosmological constant A by i?^ = 
— (d — 2) (d — 1)/2A. The parameter tq is proportional to the 
mass M of the black hole: M = {d - 2)Ad_2»'o"'Vl67r, 
where yld_2 = 27r('*-i)/Vr[(d - l)/2]. The well-known 
Schwarzschild geometry corresponds to i? — > oo. The black- 
hole horizon radius r = r_|- is the (unique) positive real root 
of/(r) =0. 

In particular, we want to focus here on the large black hole 
regime, r+/R oo. In this case, the above geometry goes 
over to a d— dimensional plane-symmetric spacetime (black 
brane) 1 16], which is also an exact solution of Einstein's equa- 
tions 1,22, 23„ .24„ .25il . The black brane has the following line 
element 

ds^ = -f(r)dt^ + y dx'dx., + —— dr^ , (2.2) 

/(^) 

where /(r) = r'^/R'^ - r'^^^ / {R^r'^-'^). The Hawking tem- 
perature of the black hole ( 12.1b is 



T = 



{d - l)r\ + [d ~ ■i)R^ 
Anr+R^ 



(2.3) 



and it reduces to T = (rf — \)r+/ AttR^ for large black holes, 
which is the Hawking temperature of the black brane ( 12.2b . 



B. Gravitational perturbations 



il perturbations in these back grou nds were con- 
IHIIIIII for d = 4 and in llills 
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sidered in Ref. |[l7t|26|,|27|] for d = 4 and in ||28l|29tl for higher 
dimensions. In a generic number of dimensions, the gravita- 
tional perturbations can be divided in three different types, the 
tensor-, vector- and scalar-type perturbations. These can all be 
reduced to a master wave equation of the form 



,d^^ 



,d^ 



where the potential V{r) depends on both the type of pertur- 
bation and on basis functions used to separate the coordinates 
of the {d — 2)-dimensional hypersurface of constant r and t. 



For instance, for tensor-type perturbations in the black-hole 
background ( 12. lb , 

V lil + d-3) (d-2)(d~4) {d-2)f' 

f - j,2 + 4^2 / + 2r ■ ^ ' 

Here, the angular number I is related to the eigenvalue of the 
hyper-spherical functions used to factor out the dependence 
on the {d — 2)-spherical angles. 

For the black brane background ( 12.21 ). the separation of triv- 
ial dimensions is achieved through the Ansatz e*'''^ and one 
ends up with the large r+ /i? limit of (12.5b . 
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The explicit form of V{r) for the vector- and scalar-type per- 
turbations can be found in Refs. |[l7[|28ll29ll . The potential for 
tensor-type gravitational perturbations is equal to the potential 
for scalar-field perturbations, so our results are also valid for 
spin-0 fields. 

We will be interested in the large-?,;; limit of the QNM 
frequencies, i.e., characteristic cj's for which the solutions to 
( 12.4b satisfy the appropriate boundary conditions. Note that 
in this limit one can formally identify q with I in the black 
hole/black brane spacetimes. Hereafter we will always refer 
to q, with the understanding that the replacement q I de- 
scribes large black holes. 

An important characteristic of classical field evolutions on 
asymptotically AdS spacetimes is the variety of choices for 
the boundary conditions at spatial infinity. In general, these 
can be Dirichlet, Neumann or Robin boundary conditions. So 
we need to establish an objective criterion before to choose a 
specific condition. In the AdS/CFT context, a natural criterion 
is such that the QNM frequencies of a certain field correspond 
to poles of two-point correlation functions of the dual operator 
in the boundary field theory li^ 30, 31, 3Z1 . 

When we consider a variable 5" such that the master equa- 
tion for the gravitational perturbations is of the form ( 12.41 ), 
in general the Dirichlet condition at spatial infinity is the 
'correct' boundary condition. There is only one exception: 
for scalar-type perturbations in four spacetime dimensions, 
the boundary condition that leads to QNM frequencies corre- 
sponding to poles of retarded correlation functions is of Robin 
type (see Refs. ll33l [34ll for a discussion in d = 4, 5 and ll36ll 
for an arbitrary number of dimensions.) 



C. Long-lived modes in the eikonal limit: analytical prediction 

In asymptotically AdS spacetimes the eikonal limit is espe- 
cially interesting, since large-q modes can be very long-lived 
lITii [l6ll . A WKB analysis suggests that for the tensor-type 
gravitational perturbations (and therefore also scalar fields) 
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and r+/R^ 1 fl%, the following asymptotic behavior holds 

2d-2 

g + n„(-±j q-^ , (2.7) 
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as g ^ oo. So large-g modes are very long-lived, and they 
could play a prominent role in the BH's response to generic 
perturbations. This is at variance with the asymptotically flat 
case, where the damping timescale is roughly constant as q 
varies. Notice also that the scaling with the BH size differs 
from that of the weakly- and highly-damped modes. 



brane. The ringdown is characterized by the decay timescale 
and ringing frequency, which can directly be extracted from 
the slope and frequency of the signal above. Equivalently, 
we characterize the QNMs by a complex frequency oj = 
ur — iuJi. Time evolution methods cannot determine very 
accurately which overtone is dominating the response, though 
experience has shown that the fundamental mode seems to be 
more excited than all others, and by definition decays much 
slower iQl]. Most importantly, the scattering of wavepackets 
shows that weakly-damped modes are excited in physically 
interesting situations. Although not the direct focus of this 
work, our numerical results show no sign of apower-law tail 
at late stages, confirming earlier predictions ifmlssll . 



B. Numerical Results 



III. QUASINORMAL FREQUENCIES 
A. Methods 




FIG. 1: Typical time-domain evolution of a gaussian wavepacket 
resulting in ringdown (top to bottom d = 4, 5, 6). Here we take 
V'(m, -30) = 0, ip{Q, v) = exp [-{v - 25)^/18] and g = 5, where 
u, V are standard null coordinates |[l9ll20ll . The signal is measured at 
r«/i? = — 1, where dr/dr, = f{r) and spatial infinity is at r« = 0. 



We use two conceptually different methods to determine the 
gravitational QNM frequencies of the spacetimes ( 12. It and 
(12.2b . and both methods yield consistent results, within the 
expected error bars. The first consists on a series expansion 
method iQJ, [T6ll . which reduces the problem to finding roots 
of a polynomial. This method is well suited to large black 
holes and black branes, though the convergence properties 
worsen for large wavenumber q. In fact, for larger dimensions 
{d > 6) and higher overtones (n > 1 or 2) the problems of 
convergence of the series solution arise even for intermediate 
wavenumber values {Rq/r+ ^ 10). 

The second method employed in this work consists on a 
direct time-evolution in these backgrounds lfl9i I20I1 . A par- 
ticular example is shown in Figure [T] which shows the time- 
development of a Gaussian wavepacket in a ?-+ /i? = 1 black 
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FIG. 2: Numerical results for the fundamental scalar-field (tensor- 
type) QNM frequencies of a r+/i? = 1 black brane. Upper panel: 
real component lor (top to bottom are d = 6, 5, 4). Lower panel: 
imaginary component uii (top to bottom for g > 10 corresponds to 
d = 4, 5, 6). Dotted lines are the analytical prediction J2.7t . cor- 
rected by a prefactor a shown in Table HH 

In Figure |2] we show numerical results for scalar-field 
(tensor-type gravitational) perturbations of a r+/i? = 1 black 
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brane. Similar results hold for large black holes. Low-g re- 
sults for d = 4, 5 are well known in the literature (See Ref. ^ 
and references therein), while the higher-dimensional cases 
are discussed in details in Ref. Isill : both ujb, and uii are 
almost independent on q (or I for large black holes) in this 
regime. For wavenumbers q S> r+/R, the qualitative be- 
havior changes: ujb_ grows linearly while loi decreases with 
q. Furthermore, it is also clear from Fig. |2]that both the sub- 
leading term in luji and the leading term in loi scale as a power 
of q (or 1). This power can be directly read from the slope of 
the curves of Figure |2] at large q. To investigate this further, 
we parameterize the numerical results by 

RujR = q + anq-'^'' , Rloj ^ aj q-f^' , (3.1) 

and we extract / by least-squares. We obtain the values 
listed in Table U where we also show the prediction by FL, 
i.e.,/3K,/ = (d-3)/(d+l). 



TABLE I: Best fit to exponents Pnj, and PL's prediction, ^qrf- 



Type 


n 


4d 




5d 




6d 






Pr Pi 


FL 


Pr Pi 


FL 


Pr Pi FL 







0.20 0.20 


0.20 


0.33 0.33 


0.33 


0.43 0.43 0.43 




1 


0.21 0.20 


0.20 


0.34 0.33 


0.33 


0.43 0.43 0.43 


tensor 


2 


0.21 0.20 


0.20 


0.34 0.33 


0.33 


0.43 0.43 0.43 




3 


0.22 0.19 


0.20 


0.34 0.33 


0.33 


0.43 0.42 0.43 




4 


0.24 0.20 


0.20 


0.34 0.32 


0.33 


0.43 0.41 0.43 







0.19 0.20 


0.20 


0.33 0.34 


0.33 


0.42 0.44 0.43 


vector 


1 


0.20 0.20 


0.20 


0.33 0.33 


0.33 


0.43 0.43 0.43 




2 


0.20 0.20 


0.20 




0.33 


- - 0.43 







0.17 0.21 


0.20 


0.30 0.31 


0.33 


0.39 0.41 0.43 


scalar 


1 


0.16 0.22 0.20 


0.31 0.35 


0.33 


- - 0.43 




2 


0.19 0.21 


0.20 


0.31 0.34 0.33 


- 0.43 



Our numerical results are consistent with a 
dependence of the characteristic frequencies, not only for the 
dimensions shown in Table|I] but also for d = 7,8 and 9. Fur- 
thermore, we computed the same modes for a 7-+/i? = 100 
black hole, and to numerical accuracy we get the same results 

2d-2 

after a rescaling by (^) ''^^ is performed. Thus, our results 
are also highly consistent with the functional dependence on 
r+, g as given by equation ( 12.7b . This agreement is nicely il- 
lustrated in Fig. |2] In a log-log plot, the analytical result pre- 
dicts a line with slope — fqrf , which overlaps very well with 
the numerical results for large q. 

We now assume the power-law behavior (I2.7l i in q and r+, 
and fit the numerical results to the following function 

2d-2 

Rlo^^"' = l+(ai^Re[n„] +m,Im[n„]) (^) q-^ , 

(3.2) 

thereby testing the prefactor in ( 12. 7b . If Eg. (12. 7b captures cor- 
rectly all of the features of these modes, then a/j w a/ w 1. 
Table summarizes our main results, with the correction 



factors to (12.7b for each dimension number d and gravita- 
tional sector The results in Table are strong indicators 
that Eg. ( 12. 71 ) does not account for the correct quantitative 
behavior of these weakly-damped modes. The results are 

TABLE II: Correction factors to the analytical formula l i2.7t . 



Type 


n 




4d 




5d 


6d 






aR 


ai 


a.R 


ai 


clr ai 




n 
U 


1.83 


1.82 




1 QQ 


4. JO 4.JJ 




1 


1.86 


1.85 


3.12 


3.10 


4.82 4.80 


tensor 


2 


1.87 


1.86 


3.15 


3.13 


4.86 4.90 




3 


1.88 


1.86 


3.17 


3.14 


4.94 4.89 




4 


1.86 


1.85 


3.18 


3.14 


4.96 4.89 







1.00 


1.02 


1.80 


1.82 


2.92 2.95 


vector 


1 


1.38 


1.38 


2.34 


2.35 


3.68 3.69 




2 


1.53 


1.53 













0.29 0.30 


0.77 


0.80 


1.59 1.64 


scalar 


1 


0.92 


0.93 


1.59 


1.65 






2 


1.20 


1.20 


2.02 


2.04 





consistent with real, overtone-independent, but dimension- 
dependent correction factors for scalar-field (tensor-type grav- 
itational) perturbations. This correction factor grows with d 
and might become dominant at large d. 

Equation ( 12.7b is not supposed to hold for vector-type and 
scalar-type gravitational perturbations, but we find it captures 
the essential qualitative behavior with r+, q. It can describe 
quantitatively the numerical results if multiplied by a real con- 
stant, which depends on the overtone ?i and the spacetime di- 
mension d. This clearly suggests a new form for n,i. 

IV. CONCLUSIONS AND OUTLOOK 

Our numerical results lend strong support to FL's predic- 
tion for the existence of long-Uved modes in the eikonal limit. 
Furthermore, the functional dependence of these modes on the 
horizon radius and momentum q is consistent with FL, but we 
show that the correct quantitative behavior is not. In partic- 
ular, if we correct their prediction by some (real) correction 
factors, listed in Table HH one can account extremely well for 
the numerical results. Taken together, our results suggest that 
large-q tensor-type (or scalar-field) quasinormal frequencies 
of black holes and black branes are described by 

2d-2 

i?w = q + an„(^) g-^, (4.1) 

where the value of a depends on d, but it is independent on 
the overtone number n. 

We also provide correction factors that would make predic- 
tion ( 12.7b describe well other types of gravitational quasinor- 
mal modes (vector and scalar). For such perturbations, the 
real constant a depends not only on d but also on n, sug- 
gesting a completely different form to n„. In any case, there 
is a simple and universal dependence on r+ and q, I in this 
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eikonal regime. Perhaps a simple interpretation in terms of 
geodesies can be given, as is done in asymptotically flat space- 
times IiI IjtIi . Clearly, more analytical and numerical studies 
are necessary to have a clear picture of the eikonal, weakly- 
damped regime of quasinormal modes of large black holes 
and black branes. A particularly interesting direction is to as- 
sess the degree to which these modes can be excited, which is 
tantamount to computing the residue of the Green function at 
the QNM pole. This is an important research topic in asymp- 
totically flat spacetime, where it allows to predict how astro- 
physical black holes respond to external sources Issll . and has 
also recently started to be explored in the gauge/gravity dual- 
ity scenario ||39i|40|l . 

Note added in proof: We have recently been informed ll4lll 
of a mistake in one integral in FL, which introduces the cor- 



rection factor a = (l/2)(d-l)(2''-2)/(d+i)inEq. This 
correction factor is consistent with all our numerical results 
for scalar fields or tensor-type gravitational perturbations. 
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